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ABSTRACT 

In  this  paper  we  discuss  the  wave-resolution  properties  of  the  Fourier  approximations  of 
a  wave  function  with  discontinuities.  It  is  well  known  that  a  minimum  of  two  points  per 
wave  is  needed  to  resolve  a  periodic  wave  function  using  Fourier  expansions.  For  Chebyshev 
approximations  of  a  wave  function,  a  minimum  of  ir  points  per  wave  is  needed  [3].  Here  we 
obtain  an  estimate  for  the  minimum  number  of  points  per  wave  to  resolve  a  discontinuous 
wave  based  on  its  Fourier  coefficients. 

In  our  recent  work  on  overcoming  the  Gibbs  phenomenon,  we  have  shown  that  the  Fourier 
coefficients  of  a  discontinuous  function  contain  enough  information  to  reconstruct  with  expo¬ 
nential  accuracy  the  coefficients  of  a  rapidly  converging  Gegenbauer  expansion.  We  therefore 
study  the  resolution  properties  of  a  Gegenbauer  expansion  where  both  the  number  of  terms 
and  the  order  increase. 
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1  Introduction 

In  [5]  (see  also  [3]),  the  issue  of  the  number  of  points  required  to  resolve  a  wave  has  been 
considered  as  a  measure  of  the  accuracy  of  a  given  scheme.  In  particular  it  has  been  shown 
that  only  two  points  per  wave  are  required  if  the  wave  function  is  approximated  by  a  Fourier 
series.  This  result  is  almost  trivial  since  the  expansion  basis  contains  the  wave  function.  It 
is  also  shown  in  [3]  that  if  the  wave  function  f{x)  =  is  approximated  by  a  truncated 
Chebyshev  expansion  then  a  minimum  number  of  ir  points  (or,  alternatively,  number  of  terms 
retained  in  the  expansion)  per  wave  is  required.  The  proof  is  based  on  the  observation  that 
if  the  number  of  terms  in  the  expansion  exceeds  this  minimum  of  tt  per  wave  then  the  error 
decays  exponentially  (see  [3,  page  35]). 

The  situation  changes  when  the  wave  function  is  discontinuous.  As  a  generic  example, 
consider  the  above  mentioned  function  f{x)  =  c*'"*  defined  on  [—1, 1]  where  uf  is  not  an  in¬ 
teger.  In  this  case,  the  truncated  Fourier  expansion  does  not  converge  at  all  in  the  maximum 
norm.  This  is  known  as  the  Gibbs  phenomenon  and  there  is  no  meaning  to  the  question  of 
number  of  points  required  to  resolve  such  a  wave. 

Recently,  however,  the  Gibbs  phenomenon  has  been  resolved.  In  [4]  it  has  been  shown 
that  the  first  N  Fourier  coefficients  of  an  analytic,  but  nonperiodic  function  contain  enough 
information  on  the  solution  to  construct  an  exponentially  convergent  Gegenbauer  series. 

Consequently  we  are  concerned  with  resolution  properties  of  the  expansion  in  Gegenbauer 
polynomials  C^{x)  to  approximate  a  non-periodic  wave  function  f{x)  —  where  u;  is  not 
an  integer.  We  first  consider  the  question  of  minimum  number  of  points  per  wave  needed 
to  resolve  this  function.  Two  separate  situations  are  considered.  The  first  involves  a  fixed 
A  in  the  Gegenbauer  expansion.  Note  that  the  Chebyshev  expansion  is  a  special  cMe  with 
A  =  0.  For  this  situation  our  results  show  that  the  minimum  number  of  points  per  wave  is 
TT,  the  same  as  in  the  special  Chebyshev  case.  In  the  second  situation  we  assume  that  7  =  ^ 
is  a  constant  in  the  Gegenbauer  expansion.  We  show  that  the  minimum  number  of  points 
per  wave  increases  with  7  at  the  boundary  a;  =  ±1  (i.e.  at  the  discontinuity  of  the  wave) 
and  decreases  with  7  at  the  center  i  =  0.  This  second  situation  corresponds  to  the  case  in 
which  the  first  N  Fourier  coefficients  of  the  function  are  known,  since  in  this  case  we  can 
reconstruct  the  Gegenbauer  expansion  for  7  =  ^  constant. 

In  Section  2  we  quote  some  results  about  Gegenbauer  polynomials.  In  Section  3  we  study 
the  resolution  properties  of  the  Gegenbauer  expansion  of  a  wave  function,  first  for  A  fixed 
and  then  for  7  =  ^  constant.  Section  4  is  devoted  to  the  discussion  of  how  many  Gegenbauer 
coefficients  can  be  obtained  accurately  from  the  first  N  Fourier  coefficients.  Here  we  have 
the  minimum  number  of  points  per  (discontinuous)  wave  to  get  exponential  decay  of  the 
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error.  Section  5  brings  numerical  confirmation  for  the  above  results. 

We  will  use  A  to  denote  a  generic  constant  independent  of  all  the  growing  parameters 
throughout  this  paper.  The  actual  value  of  A  may  be  different  at  different  locations. 


2  Preliminaries 

In  this  section  we  will  collect  some  definitions,  equalities  and  inequalities  which  will  be  used 
in  later  sections. 

We  start  by  defining  the  Gegenbauer  polynomials  G^Ca:). 

Definition  2.1  For  A  >  0,  the  Gegenbauer  polynomial  C^(x)  is  the  polynomial  of  degree  n 
that  satisfies 


y  (1  -  x^)‘‘-‘wi(x)Cji(x)dx  =0  k^n 


(2.1) 


with 


C„"(l) 


r(n  +  2A) 


(2.2) 


n!r(2A)  ' 

The  Gegenbauer  polynomial  C^(x)  achieves  its  maximum  at  the  boundary  i  =  1  [2,  page 


206] 


|C„'(I)|  <  C„'(l)  (2.3) 

and  its  value  at  the  center  x  =  0  is  given  by  [1,  page  777] 

CLm  =  (-!)"  (0)  =  0.  (2.4) 

The  Gegenbauer  polynomials  thus  defined  are  not  orthonormal.  The  norm  of  C',J(x)  is 
given  by  [2,  page  174] 


T(A)(n  +  A)' 

We  will  also  need  the  following  identity,  which  can  be  found  in  [2,  page  213] 


(2.5) 


i  y(l-i“)*-ie''“'C,'(i)(ii  =  r{A)(|^)  i'(H-A)J,+j(Ti.>).  (2.6) 

Throughout  this  paper  we  repeatedly  use  the  Stirling’s  formula:  for  any  x  >  1, 

r(x  +  l)  <  (27r)»x*+»e-‘'e»  (2.7) 

r(x-l-l)  >  (2?r)*x^+»c-*  .  (2.8) 
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We  will  need  the  following  estimates  for  the  Bessel  function  Jn{nz)  [1,  page  362] 


|7„(nz)| 


<( . 


0  <  z  <  1. 


The  following  Lemma  will  be  used  in  conjunction  with  the  estimate  in  (2.9). 

Lemma  2.1 

If  q{z)  and  p{z)  are  defined  by 


9(^) 


p{z)  = 


z'v 


where  7  >  0,  and  the  constant  0(7)  is  defined  by 


then 


c(7)  = 


n/TT^ 
1+7  ’ 


(1)  For  0  <  z  <  1,  q(2)  is  a  strictly  increasing  function  and  ^(z)  <  1; 

(2)  For  0  <  z  <  0(7),  p(z)  is  a  strictly  increasing  function; 

(3)  For  0(7)  <  z  <  1,  p{z)  a  strictly  decreasing  function  and  p(z)  >  1. 


(2.9) 


(2.10) 


(2.11) 


□ 


3  Wave  Resolution  Properties  of  Gegenbauer  Expan¬ 
sions 

Consider  the  non-periodic  wave  function 


f{x)  =  e'^^  (3.1) 

where  a;  >  0  is  not  an  integer.  The  Gegenbauer  coefficients  of  this  function  are  defined  by 

/^(O  =  ^  -  x*)^-ie’^C',^(x)dx  (3.2) 

with  hf  given  by  (2.5). 

Our  objective  is  to  find  a  lower  bound  on  the  ratio 


where  u;  is  the  number  of  waves  in 
expansion 


_ 


(3.3) 


(3.1)  and  m  is  the  number  of  terms  in  the  Gegenbauer 
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(3.4) 


1=0 

such  that  the  approximation  error  using  the  expansion  (3.4)  is  exponentially  small  when 
m  — »  c».  The  ratio  r  defined  in  (3.3)  is  usually  called  the  number  of  points  per  wave.  It  is 
literally  the  number  of  points  per  wave  for  the  collocation  case  and  can  be  called  number  of 
modes  per  wave  for  the  Galerkin  expansion  (3.4)-(3.2).  We  can  then  define  the  regularization 
error  RE(A,m,r,  x)  to  be 


RE(A,  m, r, x)  =  |/(x)  -  /^(z)|  (3.5) 

and  ask  the  question  of  finding  a  lower  bound  Tq  =  ro(A,  x)  such  that  the  regularization  error 
(3.5)  is  exponentially  small  for  r  >  tq  when  m  — ^  oo. 

According  to  the  identity  (2.6)  and  the  definitions  (3.2)-(3.4),  the  regularization  error 
(3.5)  can  be  explicitly  expressed  as 

flE(A,m,r,x)=  r(A)  f— f]  i'{/ +  A)J,+,  (— )  C,V)  •  (3.6) 

first  we  have  the  following  theorem  for  the  case  with  fixed  A. 

Theorem  3.1  If  A  is  fixed,  then  ro  =  is  a  lower  bound  for  the  number  of  points  per  wave 
to  obtain  exponentially  small  regularization  error  for  all  —1  <  x  <  1. 

Proof:  Assume  r  >  tt.  Denote  zj  =  and  z  =  For  I  >  m  we  have  z/  <  z  <  1  and 

hence  q(zi)  <  q(z)  <  1  according  to  Lemma  2.1.  Take  s  =  j  and  m  big  enough  so 

>171.  If  we  define 

B(()  =  ((+A)(,(x))'+^iy.^  (3.7) 

we  have,  for  /  >  m 


B{1) 


(/+H-A)(/  +  2A) 

(/  +  A)(/  +  1) 


q{z)  <  s  <  1 . 


(3.8) 


We  can  then  start  from  the  explicit  formula  (3.6)  and  absorb  all  the  A  dependent  terms 
into  the  generic  constant  A  to  obtain 


RE{X,m,r,x)  < 
< 


Am-^  E  {I  AX)  Jl+J—)  C,')!) 

I=m+I  V  r  / 

l=m+i  *• 
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< 

Am  ^  5^  . 

< 

1  —  s 

< 

{q(z))”' 

where  we  have  used  (2.3)  in  the  first  inequality  (2.9)  and  (2.2)  in  the  second  inequality,  the 
definition  of  B{1)  in  (3.7),  the  fact  that  q{zi)  <  q{z)  in  the  third  inequality,  and  (3.8)  in  the 
fourth  inequality.  This  finishes  the  proof  because  q{z)  <1.  □ 


The  result  of  Theorem  3. 1  is  not  surprising;  for  A  =  0  this  is  simply  the  known  result  for 
the  Chebyshev  case  [3]. 

For  fixed  A,  the  estimates  are  essentially  the  same  for  the  boundary  point  a:  =  1  or  for 
the  center  point  1  =  0,  since  there  is  only  an  algebraic  difference  between  (^/^(l)  and  ^^(0). 
However,  this  algebraic  difference  becomes  bigger  when  A  increases. 

For  our  purpose  we  are  more  interested  in  the  case  A  ~  m.  This  is  because  we  are 
interested  in  the  situation  where  the  finite  Fourier  series  of  f{x)  is  given.  In  [4]  we  proved 
that  it  is  possible  to  recover  uniform  exponential  accuracy  from  the  finite  Fourier  series  of  a 
non-periodic  analytic  function  through  the  use  of  Gegenbauer  expansions  with  A  ~  m. 

Let  us  now  assume  that  7  =  ^  is  a  constant.  We  have  then  the  following  theorem. 
Theorem  3.2  If  =  —  is  a  constant,  then 

(I)  The  regularization  error  /2£J(7m,m, r,  1),  at  the  boundary  x  =  1,  is  exponentially 
small  if  r  >  max(ri,r2)  where  rj  is  the  unique  solution  of 


(n(l  +7)) 


(2e7(l+7))^^^^(7) 


,n(l+7)/  (1+27)^+2t' 

in  the  region  rj  >  (3-9)  ^  solution  (i.e.,  if  p(c(7))  >  ^1(7)),  or  n 

otherwise.  Also,  rj  is  the  unique  solution  of 


(3.9) 


yl+J-r 


('•1(1  +  7))  1+27’ 


,-■1(1+7)/  ^ 

(2)  The  regularization  error  RE{-ym,m,r,0),  at  the  center  x  =  0,  is  exponentially  small 
if  r  >  max(r3,r4)  where  rs  is  the  unique  solution  of 


P 


( 


TT 

r3(l  +  7) 


) 


(e(l  -H7))^ 

(Id- 27)*+"' 


Q3{j) 


(3.11) 


within  the  region  ra  >  if  (3.11)  has  a  solution  (i.e.,  if  p(c(7))  >  (^3(7)),  or  ra  = 

otherwise.  Here  r^  is  the  unique  solution  of 
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''(r4(l+7))“sAT^' 

The  functions  p{z)  and  q{z)  in  (3.9)-(3.12)  are  defined  by  (2.10)  in  Lemma  2.1  and  the 
constant  0(7)  is  defined  by  (2.11). 

Proof:  Denote  zi  =  and  z  =  ^  we  have  zi  <  z  <  0(7)  when 

^  ^  hence  q{zt)  <  q{z)  <  according  to  Lemma  2.1  and  the  condition  (3.10). 

Take  s  =  which  yields  (1  +  2'Y)q(z)  <  s  <  1.  Hence  we  can  take  m  big  enough 

so  that  for  I  >  m, 


{1+1  +  ym){l  +  27m) 

(/  +  7m)(/+l) 


-  ('  +  rrW) 


(1  +  2'y)q{z)  <  s  <  1 . 


If  we  now  define 


we  have  for  /  >  m 


^(0  =  (^  +  7m)(g(z)) 


B{1  +  1)  _  (/  +  1  +  7m)(/  +  27m) 


l+ym  r(/  +  27m) 


m  (/  +  7m)(/+l) 

We  can  then  start  from  the  explicit  formula  (3.6)  to  obtain 


q{z)  <  s  <  1 


(3.13) 


(3.14) 


RE{'im,m,r,\)  <  r(7m)  +  7^)  J,+-^  C7^(l) 

^  r(7m)  /2r \^B(m) 

~  r(27m)  \7rm/  1  —  5 


where  we  have  used  (2.9)  and  (2.2)  in  the  second  inequality,  the  definition  (3.13)  and  the 
fact  q(zi)  <  q(z)  in  the  third  inequality,  (3.14)  in  the  fourth  inequality  and  Stirling’s  formula 
(2.7)-(2.8)  in  the  last  inequality. 
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According  to  condition  (3.9)  and  Lemma  2.1,  the  number  is  strictly 

less  than  one  for  r  >  rj.  This  finishes  the  proof  for  x  =  1. 

The  proof  for  the  center  point  x  =  0  is  similar  and  is  thus  omitted.  □ 


Figure  1:  Minimum  number  of  points  per  wave  as  a  function  of  7.  Left:  minimum  r  =  ^ 
to  obtain  exponentially  small  regularization  error;  Right:  minimum  R  =  ^  to  obtain  expo¬ 
nentially  small  resolution  error  from  Fourier  partial  sum. 

In  Figure  1  (left)  we  show  the  curves  of  number  of  points  per  wave,  r  =  ^,  versus  7, 
both  for  resolving  the  boundary  x  =  1  (upper  line)  and  for  resolving  the  center  x  =  0  (lower 
line).  These  were  obtained  from  the  results  (3.9)-(3.12)  in  Theorem  3.2.  We  can  see  that 
the  top  curve  for  x  =  1  is  an  increasing  function  of  7,  while  the  bottom  curve  for  x  =  0  is  a 
decreasing  function  of  7  for  7  >  0.37.  At  7  =  1,  we  need  approximately  7.03  points  per  wave 
to  resolve  the  boundary  and  3.51  points  per  wave  to  resolve  the  center.  At  7  =  5,  this  two 
numbers  change  to  8.11  and  2.33,  respectively.  Figure  1  (right)  is  related  to  the  truncation 
error  to  be  discussed  in  next  section. 


4  The  Truncation  Error  from  Finite  Fourier  Series 

In  this  section  we  consider  the  situation  that  /(x)  =  is  not  known  and  only  the  first 
2N  +  1  Fourier  coefficients  f{k]  are  given.  We  would  like  to  recover,  within  exponential 
accuracy  at  all  points  in  —1  <  x  <  1,  the  finite  Gegenbauer  expansion  /^(x)  of  /(x)  as 
defined  by  (3.4).  This,  together  with  the  results  in  the  previous  section,  will  establish  the 
number  of  points  per  wave  for  resolving  discontinuous  waves  based  on  its  Fourier  partial  sum 
through  Gegenbauer  polynomials.  The  result  in  this  section  parallels  the  general  theorems 
in  [4]. 
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The  Fourier  partial  sum  of  the  discontinuous  wave  function  f{x)  =  is  given  by 


/«(!)  =  E  hty 

k=-N 

where  the  Fourier  coefficients  f{k)  are  defined  by 


Assume  that  the  first  2N  +  1  Fourier  coefficients  f(k),  |A;|  <  N  and  hence  also  the 
Fourier  partial  sum  (4.1)  are  known  but  the  function  f{x)  is  not.  We  thus  do  not  know  the 
exact  Gegenbauer  coefficients  f^{l)  of  f(x)  defined  by  (3.2),  but  only  the  approximate  ones 
obtained  from  the  Fourier  partial  sum  /yv(^)»  which  we  denote  by  g^{l) 

1  f  y.  _ A...  /  \  Vyv.  /yy.\ 


9\l)  =  ^  Jjl-x^y-^Mx)Ct{x)da 


Notice  that  g^(l)  depends  on  N. 

We  now  define  the  truncation  error  TE{\,m,  N,x)  to  be 

m 

TE(X,m,N,x)=  E(/'(0-j‘(0)C,"(l)  ■  (4-4) 

1=0 

Notice  that  u;  or  r  =  ^  is  not  an  explicit  parameter  in  our  definition  of  the  truncation 
error.  The  estimates  we  obtain  later  will  be  uniformly  valid  for  all  uj. 

In  the  next  theorem  we  shall  bound  the  truncation  error  in  terms  of  N,  the  number  of 
given  Fourier  coefficients,  m  the  number  of  terms  in  the  Gegenbauer  expansion,  and  A. 

Theorem  4.1  The  truncation  error  at  the  boundary  x  =  1  satisfies  the  estimate 


TE{X,m,N,\)  <  A 


(m  -f  A)r(m  +  2A)r(A)  /  2 


(m-l)!r(2A) 

At  the  center  x  =  0,  the  truncation  error  satisfies 


(-)  ■ 
\xNJ 


TE(X,m,N,(l)  <  A 


(m  +  A)r(f  +  A)  /  2 


r{?) 


(-)  ■ 
\wnJ 


If  A  =  ym  and  m  =  where  7  and  ^  are  positive  constants,  then  the  truncation  error 
satisfies 


TE{y0N,0N,N,\)  <  AN^ 


,  (27re7)T'  ) 


at  the  boundary  x  =  1  and 


at  the  center  a:  =  0. 
Proof;  Since 


TE{^0N,ffN,N,O)  <  AN^  j 


m  -  Mx)  =  ^  i(ny' 

\n\>N 


0N 


(4.8) 


(4.9) 


and  |/(n)|  <  1  according  to  the  definition  (4.2)  with  f{x)  =  we  can  estimate  the 

truncation  error  as  follows 


=  (m  +  1)  max 

0</<m 


TEiX,m,N,x)  <  (m+1)  max  |(/^(()-S*(0)C,\x)| 

Y:  nn)Ct(x)  (i  /’  (1  -  xy-ie‘'"''CKy)dv) 

S  /(n)C,\x)r(A)  (— )'  i'(/  +  X)J,Mxn) 

|n|>\ 

<  Amr(A)(^)  max  (I  +  A)  lC‘(x)|  . 


=  (m  +  1)  max 

0</<m 


(4.10) 


In  the  second  step  we  have  used  (4.9),  the  definition  of  f^{l)  in  (3.2)  and  that  of  g^{l)  in 
(4.3),  in  the  third  step  we  have  used  the  equality  (2.6),  and  in  the  last  step  we  have  used 
the  facts  |/(n)|  <  1  and  |J,/(x)|  <  1  for  all  x  and  >  0  [1,  page  362]. 

For  the  boundary  point  i  =  1  we  can  then  proceed  as  follows 

TE{X,m,N,\)  <  x4mr(A)^^^  max  (/ +  A)C,'^(1) 


=  i4m 


r(A) 


( 


A-l 


max 


(/  +  A)r(/  +  2A) 


<  A 


r(2A)  \irNJ  0<Km  /! 
(m  +  A)r(m  +  2A)r(A)  /  2  ' 


(m-l)!r(2A)  \irNJ 

where  in  the  second  step  we  have  used  the  formula  (2.2)  for  C^{1),  and  in  the  last  step 
we  have  used  the  fact  that  is  an  increasing  function  of  /.  The  result  in  (4.5)  is 

thus  proven.  Some  simple  algebra  and  the  use  of  Stirling’s  formula  (2.7)-(2.8)  easily  produce 
(4.7).  The  estimates  (4.8)  and  (4.10)  at  the  center  point  i  =  0  can  be  obtained  in  a  similarly 
fashion  and  the  detail  is  thus  omitted.  □ 

Since  our  estimates  of  the  truncation  error,  (4.7)  and  (4.8),  do  not  depend  on  the  wave 
number  u;  or  r  =  — ,  the  minimum  number  of  points  per  wave  for  exponential  convergence. 
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defined  by  is  obtained  by  using  the  ^  derived  in  Theorem  3.2  and  the  largest 

0  in  (4.7)  or  (4.8)  such  that  the  factors  on  the  right  hand  sides  which  are  being  raised  to 
power  N  are  still  less  than  one.  This  gives  the  following  theorem. 

Theorem  4.2  Assume  that  the  Fourier  coefficients 

m  =  5  /_'  /(x)e--‘'Vx 

of  the  function  f{x)  =  are  known  for  —N  <  k  <  N.  Let  0  <  /  <  m  be  the 

Gegenbauer  expansion  coefficients  of  /Ar(i)  =  EtL-N given  by  (4.3).  Define  the 
number  of  points  per  wave  by 


U} 

and  the  resolution  error  by 


(4.13) 

Then,  if  A  =  7m  and  m  =  l3N^  we  have  the  following  results. 

( 1 )  At  the  boundary  x  =  1 ,  if  one  uses 


E(A,m,yV,  jR,  x)  = 


/(i)  - 


1=0 


(4.12) 


0< 

and  r  =  ^  >  max(ri ,  rj)  as  is  derived 
tially  small  and  bounded  as 


_  r,  /  X 

- - r-7XT7  =  '^1(7) 

1+27)^+* 


(4.14) 


in  Theorem  3.2,  then  the  resolution  error  is  exponen- 


E(-rl)N,  RN,  N,R,l)<  A  (W’i?  +  6g) 


(4.15) 


where 


V  pire^)"'  j  ’  "  V(2e7(l +7))''’^(r(l +7)/j 


<  1 


The  number  of  points  per  wave  R  =  ^  is  estimated  by 

2max(ri,r2) 


R> 


Bih) 


(4.16) 


where  T],  r2  are  defined  in  Theorem  3.2  and  ^1(7)  is  defined  by  (4.14). 
(2)  At  the  center  x  =  0,  if  one  uses 


/3< 


ve 


(1+27)*+' 


(4.17) 


and  r  =  ^  >  inax(r3,  r4)  as  is  derived  in  Theorem  3.2,  then  the  resolution  error  is  exponen¬ 
tially  small 


E(-iffN,  ffN,  N,  R,0)<  A  +  cjf) 


where 


CT=l - ^  ^  " - j  <1 

(Tre)^  / 


Cr  = 


(1+27)> 


(e(l+7)) 

The  number  of  points  per  wave  R  =  ^  is  estimated  by 

2max(r3,r4) 


?+->  /  TT  ' 

)r'’wi+7), 


R> 


^2(7) 


where  r3,  r4  are  defined  in  Theorem  3.2  and  is  defined  by  (4.17). 
Proof:  This  is  simply  a  combination  of  Theorems  3.2  and  4.1. 


(4.18) 


<  1. 


(4.19) 


□ 


In  Figure  1  (right)  we  show  the  curves  of  number  of  points  per  wave,  R  =  ~,  versus  7, 
both  for  resolving  the  boundary  a:  =  1  (upper  line)  and  for  resolving  the  center  x  =  0  (lower 
line),  obtained  from  the  results  (4.16)  and  (4.19)  in  Theorem  4.2.  We  can  see  that  the  top 
curve  for  X  =  1  achieves  its  minimum  at  around  7  =  0.9.  and  the  bottom  curve  for  x  =  0 
is  an  increasing  function  of  7.  We  emphasize  that  these  results  are  obtained  with  different 
;5  for  X  =  1  and  x  —  0  (from  (4.14)  and  (4.17)  respectively).  In  practice  a  single  ^  should 
be  used  since  one  would  like  to  resolve  both  the  boundary  and  the  center  simultaneously.  If 
the  ^  for  X  =  1,  given  by  (4.14),  is  also  used  for  x  =  0,  the  minimum  number  of  points  per 
wave,  R,  to  resolve  the  center  x  =  0,  would  be  described  by  the  middle  curve  in  Figure  1 
(right).  For  the  single  0  chosen  according  to  (4.14),  at  7  =  1,  a  minimum  of  22.2  points  per 
wave  is  needed  to  resolve  the  boundary  x  =  1,  and  a  minimum  of  11.1  points  per  wave  is 
needed  to  resolve  the  center  x  =  0. 


5  Numerical  Results  and  Conclusions 

In  this  section  we  perform  numerical  calculations  to  demonstrate  the  theory  developed  in 
previous  sections.  We  use  the  discontinuous  wave  function 

/(x)  =  cos(7ra;(x  -|- 1))  (5,1) 

for  various  wave  numbers  u  and  report  both  the  regularization  error  defined  by  (3.5)  and 
the  resolution  error  defined  by  (4.13). 
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We  implement  the  method  in  the  following  way.  The  exact  Gegenbauer  coefficients  of 
f{x),  which  are  needed  for  the  regularization  error,  are  computed  using  (2.6).  The  approx¬ 
imate  Gegenbauer  coefficients  defined  in  (4.3),  which  are  needed  for  the  resolution 

error,  are  computed  using  the  following  formula 

9‘(i)  =  W(0)  +  r(A)i'((  +  A)  2  m  (5.2) 

o<|fc|<Ar  VTrfc/ 

where  f{k)  are  the  Fourier  coefficients  of  f{x)  defined  by  (4.2).  This  formula  can  be  easily 
derived  from  the  definition  of  g^{l)  in  (4.3)  and  the  integration  formula  (2.6).  We  compute 
the  Bessel  function  J^{x)  using  an  IMSL  routine.  The  approximation  to  f{x)  is  obtained  by 
directly  summing  (3.4)  for  the  regularization  error,  or  by  directly  summing 

m 

(5.3) 

/=o 

for  the  resolution  error.  The  Gegenbauer  polynomials  Ci'{x)  are  computed  by  the  formula 

=  t  Trw*(f-~.)!r(A)* 

which  can  be  found  in  [1,  page  175]. 

We  remark  that  the  implementation  techniques  described  above  are  subject  to  roundoff 
effects  for  large  A  and  m.  A  better  way  of  implementing  the  method  through  Chebyshev 
polynomials  is  currently  under  investigation. 

In  Figure  2  we  show  the  errors,  in  a  logarithm  scale,  at  the  discontinuity  x  =  1  for  u;  =  1 .4, 
2.4,  3.4  and  4.4.  We  choose  7  =  ^  =  1  in  this  illustration.  On  the  left,  the  regularization 
error  RE{m,m,^,l)  is  shown  as  a  function  of  m,  the  number  of  terms  retained  in  the 
Gengenbauer  expansion.  On  the  right,  the  resolution  error  E{^,  y,yV,  1)  is  shown  as  a 
function  of  2iV,  to  total  the  number  of  terms  in  the  Fourier  expansion.  Here  we  take  =  0.5 
which  satisfies  (4.14).  We  can  see  that  the  errors  are  order  0(1)  until  r  =  ^  or  /?  =  ^ 
reaches  the  critical  values  obtained  in  Theorem  3.2  and  in  Theorem  4.2,  after  which  the 
errors  drop  exponentially.  Figure  3  shows  the  same  result  but  for  the  center  point  x  =  0. 
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Figure  2:  Logarithmic  error  at  the  discontinuity  a:  =  1  for  (5.1)  with  a;  =  1.4,  2.4,  3.4  and 
4.4.  Here  7  =  ^  =  1.  Left;  the  regularization  error  RE{m,ni,^,\)  as  a  function  of  m; 
Right:  the  resolution  error  E{^,  y,  N,  1)  as  a  function  of  27V,  ^  =  0.5. 


Figure  3:  Same  eis  Figure  2  but  for  the  center  point  a:  =  0. 
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